Morphology and Surface Brightness Evolution 
of z~ 1.1 Radio Galaxies 



Nathan Roche, Stephen Eales 

Department of Physics and Astronomy, University of Cardiff, 
P.O. Box 918, Cardiff CF2 3YB, Wales. 



Steve RawUngs 

Department of Astrophysics, University of Oxford, 
Nuclear and Astrophysics Laboratory, Keblc Road, Oxford 0X1 3RH, England. 



Abstract 

We use /C-band (2.1 fim) imaging to investigate the angular size and morphology of 10 6C radio 
galaxies, at redshifts 1 < z < 1.4. Two appear to be undergoing mergers, another contains two intensity 
peaks aligned with the radio jets, while the other seven appear consistent with being normal ellipticals 
in the if -band. 

Intrinsic half-light radii are estimated from the areas of each radio galaxy image above a series of 
thresholds. The 6C galaxy radii are found to be significantly smaller than those of the more radioluminous 
SCR galaxies at similar redshifts. This would indicate that the higher mean if -band luminosity of the 
SCR galaxies results from a difference in the size of the host galaxies, and not solely from a difference in 
the power of the active nuclei. 

The size-luminosity relation of the z ^ 1.1 6C galaxies indicates a 1.0-1.8 mag enhancement of the 
rest-frame i?-band surface brightness relative to either local ellipticals of the same size or FRII radio 
galaxies at z < 0.2. The SCR galaxies at z ~ 1.1 show a comparable enhancement in surface brightness. 
The mean radius of the 6C galaxies suggests that they evolve into ellipticals of L ^ L* luminosity, and is 
consistent with their low redshift counterparts being relatively small FRII galaxies, a factor ~ 25 lower 
in radio luminosity, or small FRI galaxies a factor of ~ 1000 lower in radio luminosity. Hence the 6C 
radio galaxies may undergo at least as much optical and radio evolution as the SCR galaxies. 



1. Introduction 

Powerful radio galaxies at z < 0.5 appear similar in morphology and profile to normal giant ellipticals 
(Lilly, MacLean and Longair 198S), with luminosities in the range from L* to ~ 2 mag above L* and 
no obvious correlation between radio and optical luminosity (Laing, Riley and Longair 198S; Owen and 
Laing 1989; Owen and White 1991). The relative uniformity of the optical properties of radio galaxies, 
combined with the ease of detection at higher redshifts, means that they provide a useful probe of galaxy 
evolution. Lilly and Longair (1984) measured if -band (2.2/im) magnitudes for SCR radio galaxies over 
a very wide range of redshifts (the SCR catalog is a flux limited sample containing the radio galaxies 
with apparent flux exceeding 10 Jy at 178 MHz). The K-z relation of these galaxies, which showed little 
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scatter out to z ^ 2, indicated an evolutionary brightening with redshift, consistent with a high formation 
redshift {z ~ 5) followed by passive stellar evolution. 

The radio luminosities of the most powerful radio galaxies also evolve. Laing, Riley and Longair 
(1983); using a V/Vmax test, showed that their decrease in hmiinosity with time was comparable to the 
evolution of quasars. Padovani and Urry (1992) studied the radio evolution of the 3CR catalog in more 
detail, deriving a radio luminosity function and obtaining a best fit for an exponential Lrad evolution 
with timescale t = O.lTj^Q Qg-fft^'^ i-c. t = S.B^^q ^ Gyr. It must be noted that this exponential decrease 
describes the evolution of the radio luminosity function, rather than that of individual radio galaxies. The 
galaxies might, for example, undergo short bursts of radio emission, separated by much longer periods of 
radio quiescence, with subsequent bursts decreasing with time in their peak radio intensity. 

Higher redshift radio galaxies differ in appearence from normal ellipticals, especially at shorter wave- 
lengths. Rigler et al. (1992) studied 13 SCR galaxies at 0.8 < z < 1.3 and found that they could be 
decomposed into 'active' and 'passive' components. The 'active' components were elongated, blue (ap- 
proximately /,y cx i^"), and aligned with the radio axis. Best et al. (1996), using WFPC2 imaging, found 
that in radio galaxies with a relatively small separation between the radio lobes, the blue components 
consisted of strings of several bright knots, whereas in radio galaxies with radio lobes separated by more 
than ^ 200 kpc, the blue components were much more compact. This suggested an evolutionary sequence 
in which the elongated optical components were formed as the outward passage of the radio jets triggers 
large-scale star-formation, and later fall back towards the galaxy centres. However, the aligned compo- 
nents of at least some z ~ 1 radio galaxies show some polarization (Tadhunter et al. 1992; Leyshon and 
Eales 1997), suggesting that dust-reflected light from the active nucleus is also important. 

The aligned components are less prominent at longer wavelengths, typically contributing only ~ 10% 
of the total flux at \rest ~ l^ni (i-G- the observed if -band), but many of the z ^ 1 3CR galaxies still show 
some alignment between their radio axes and iiT-band isophotes (Dunlop and Peacock 1993; Ridgway and 
Stockton 1997). Most of the red light is produced by the 'passive' components, which are redder, diffuse 
and symmetric, apparently unaligned with the radio axis, and appear to be underlying giant ellipticals. 

In the infra-red ff-band (1.65/im), Rigler and Lilly (1994) found the profile of the z = 1.18 radio 
galaxy 3C65 to be well-fitted with a pure bulge (r^) profile of effective radius = 1-77 ± 0.18 arcsec. 
Best et al. (1997), observing in the ii'-band and with the HST WFPC2, confirmed re = 1-7 arcsec 
for 3C65, and found all 21 of a sample of z 1 3CR galaxies to be well-fitted by bulge models with 
0.7 < fe < 3.9 arcsec. The size-luminosity relation of these galaxies was consistent with that of of local 
ellipticals with the ~ 1™ of brightening expected from passive evolution. Only two of the 21 galaxies 
contained substantial (> 10% of total fiux) nuclear point-sources components, although several showed 
some evidence of an excess of surface brightness above a bulge- model profile at large radii (r > 35 kpc), 
like that in cD galaxies, which might be evidence for a rich cluster environment. 

In this paper we investigate the ii'-band morphology and radii of 10 galaxies at similar redshifts but 

with more moderate radio luminosities, selected from the 6C survey which has a fiux limit about six times 
fainter than 3CR. There are some unsolved problems concerning the differences between the 3CR and 6C 
galaxies. Firstly, 6C radio galaxies at 0.6 < z < 1.8 were found to be on average 0.7 mag less luminous 
in the observed i^T-band than 3CR galaxies at the same redshifts (Eales and Rawlings 1996; Eales et al. 
1997), but it was unclear whether this was due solely to a greater near-IR fiux from the SCR nuclei, or 
to a greater mass of stars in the SCR galaxies. We shall investigate whether the differences in the radio 
and optical luminosities of SCR and 6C galaxies at z ~ 1 are related to any differences in the old stellar 
component visible in the if-band. Secondly, as the SCR and 6C galaxies at lower redshifts are similar 
in mean JT-band luminosity, the K-z relations of the SCR and 6C catalogs are different in slope, with 
the latter consistent with no evolution. Hence if high redshift 6C galaxies evolve into the low redshift 6C 
galaxies, they undergo surprisingly little change in /i'-band luminosity, but it is also possible that their 
low redshift counterparts are galaxies of lower K and radio luminosity. By determining a size- luminosity 
relation for the 6C galaxies and comparing with both z ^ 1 SCR galaxies and low-redshift radio galaxies 
(from Owen and Laing 1989), we may compare the surface brightness evolution of the two catalogs and 
determine whether the 6C galaxies are truly non-evolving. 

We assume Hq = 50 km s~'^Mpc~'^ and qo = 0.05 throughout. Section 2 describes the observational 
data. Section 3 its analysis and photometry. In Section 4 we describe the estimation of the size and 
morphological type of the galaxies, and in Section 5 discuss the appearence of each individual 6C galaxy. 
In Section 6 we present the size-luminosity relation and compare with SCR galaxies at similar redshifts 
and local radio galaxies, and in Section 7 discuss the implications for radio galaxy evolution. Section 8 
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summarizes the main conclusions. 



2. Observational Data 

The observational dataset of this paper consists of 10 images, each centred on a 6C radio galaxy with 
a known spectroscopic redshift in the 1.0 < z < 1.4 range {zmean = 1.11). This data was previously used 
by Bales et al. (1997) in a study of the i^-band luminosity evolution of radio galaxies. The 6C galaxies 
considered here have 151 MHz fluxes of 2.2-4.4 Jy, corresponding at these redshifts to L(151MHz) ~ 10^^'^ 
W Hz^^. AUington-Smith (1982) and Bales (1985) give more information on the radio source properties, 
Bales et al. (1997) and references therein detail the identification of the optical counterparts, and Rawlings 
et al. (1997) describe the spectroscopic observations. For comparison we include in our dataset an image 
of a radio-loud QSO, 5C6.8, which lies in the same redshift range as the other objects. Table 1 lists the 
radio galaxy co-ordinates and redshifts for all 11 images. 

These fields were observed over the period 18-20 January 1995 using the Redeye camera on the 
Canada- Prance Hawaii Telescope (CFHT). Redeye is a 256 x 256 HgCaTe mosaic with a pixelsize of 0.50 
arcscc, covering approximately 2x2 arcmin. Bach field was imaged nine times, with the camera being 
offset by 8 arcsec between exposures in a 3 x 3 grid. The total exposure times varied from 26 to 47 minutes 
per field. The observations were made in the the K' band, centred on a wavelength of 2.1/xm where there 
is some reduction of sky background compared to that in the standard 2.2 /im K-hand (Wainscoat and 
Cowie 1992). All 11 images had good seeing, estimated as FWHM ~ 0.9-1.2 arcsec. 

3. Data Analysis and Galaxy Detection 

After initial image processing using the IRAF package (see Bales et al. 1997), the Starlink PISA 

(Position, Intensity and Shape Analysis) package, developed by M. Irwin, P. Draper and N. Eaton, was 
used to detect and catalog the objects on each field. Object were detected on the basis that they exceed 
an intensity threshold of 1.5cr above the background noise (a being separately determined by PISA for 
each of the frames) in at least 6 connected pixels. Fluxes were estimated using the 'total magnitudes' 
option within PISA, which counts photons above the sky background level in elliptical apertures about 
the centroid of each detection, with the aperture size and shape being fitted by a curve-of-growth analysis 
to the intensity profile of each object. 

Observations of iC-band standard stars provided a photometric zero-point for the Redeye data of 
K = 21.28 for 1 count sec~^. However, due to the offset between the K and K' passbands, this zero point 
will be exact only for objects with the same K' — K as the standard stars. Radio galaxies with elliptical- 
like colours will at 1 < ^ < 1.4 be redder than the calibration stars by an estimated A{K' — K) ^ 0.13 
mag (see Bales et al. 1997), so to correct for this we add A{K) = —0.13 mag to the magnitudes, adopting 
a zero-point of K = 21.15 for 1 count sec~^. 

PISA also produces a list of 8 areas for each detected object. The first area is simply the number of 
pixels where the intensity equals or exceeds the chosen detection threshold, while sizes 2-8 are the areas 
above higher thresholds of Ithr + 2-'''"^ coimts, where Ithr is the threshold and j an integer from 2 to 8. 
For these images, Ithr ~ 50 counts so the 8th threshold would be ~ 1074 counts and ~ 7.5 mag above 
the detection threshold. 

We are only concerned here with the radio galaxies on each image; the properties (e.g. clustering) of 
the other galaxies detected will be discussed in a separate paper. Figures 1 and 2 show grcyscalc images 
and contour plots of 15 x 15 arcsec areas centred on each radio galaxy. Table 1 lists the K-band total 
magnitudes of the radio galaxies, as measured from this data (with errors estimated from the best-fitting 
models described in Section 4 below). These are consistent with the magnitudes of Bales ct al. (1997). 
Table 1 also lists the detection thresholds of the 11 images in terms of K mag arcsec"^. Two of the 
radio galaxies were detected by PISA as double objects on these images, 6C011-I-36 with components of 
apparent separations 1.7 arcsec and 6C1256-I-36 with components of separation 2.5 arcsec. These pairs arc 
still connected at the detection threshold so are only split when PISA is run with the 'deblend' option. In 
this analysis we have used PISA without deblending, so these radiogalaxies are classed as single objects, 
and the derived magnitudes and areas will be those of both components together. 

4. Estimating GalcLxy Sizes and Profiles 
4.1 Method 
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Most of the radio galaxies appear extended on these images. To estimate their intrinsic angnlar sizes 
and distinguish between disk and bulge profiles, we compare their areas above 8 thresholds as given by 
PISA with those measured for simulated galaxy profiles. We consider three types of source profile 

(i) A point source. 

(ii) An exponential disk, with surface brightness ^ oc ex^{—r /r^xp)- We consider exponential scale- 
lengths in the range 0.07 < r^xp < 1-00 arcsec. The half-light radius Vhi — LdiSr^-xp. 

(iii) A bulge (elliptical galaxy) profile, with surface brightness fx oc exp — 7.67(r/re)'''^^. We consider 
effective radii (re) in the range 0.07 < re < 3.00 arcsec. For this profile rni = Vg- 

Each model profile was generated on a 0.1 arcsec pixel grid, and then rebinned into 0.5 arcsec pixels. 
The model profiles were arranged in a grid pattern on a 'simulation image', with each model profile 
represented 12 times with slightly different positional offsets relative to the 0.5 arcsec pixel grid. The 
seeing profile was estimated for each image by fitting (using a PISA routine) a combined Gaussian- 
exponential-Lorentzian to a few objects identified as stars (of similar apparent magnitude to the radio 
galaxies). The simulation image was then convolved with the model stellar profile and normalized so that 
the model profiles had the same total intensity as measured for the radio galaxy. Gaussian noise was 
then added to the simulation image, with a equal to the sky noise as measured on the real data. 

The simulation image was then analysed using PISA, with the same detection threshold as used for 
the radio galaxy data. For each model profile, this gave 12 sets of areas above the 8 thresholds, with 
some scatter between the 12 due to the modelled sky noise and the positional offsets between the models. 
These 12 sets of areas were averaged to give a set of 8 PISA areas corresponding to each model profile, 
with error bars calculated from the scatter. For eac;li radio galaxy wt; then compared the observed set of 
8 areas with those measured for each of the models, using a test with the simulation errors. 

Table 2 shows the results of this comparison for each object, with for the point-source models, the 
exponential scalelength (rgxp) and of the disk-profile model which best fits the observed set of 8 areas, 
i.e. which gives the smallest x^, and the effective radius (re) and of the best-fitting bulge model. 

For a x^ test with 7 degrees of freedom, we would expect a 'perfect' model to on average give 
X^ — 6.4, and would estimate ilcr errors on Vexp and as the change in the model radius above or below 
the best-fitting model which increases x^ by 8.16 above its minimum value. However, in this analysis 
the 8 areas given by PISA are not entirely independent; a large noise fluctuation could change a pixel's 
value by more than one threshold interval, thus altering more than one of the areas. To take this into 
account, we also x^-test the observed areas above thresholds against each of the 12 representations of 
the best-fitting model in our simulation image, obtaining 12 values of x^- The scatter in these 12 values 
of x^ gives the ilcr errors (listed in Table 2) produced by noise on the best-fit model's x^- The errors 
on the best-fit and rexp can be estimated as the change in these radii which increases x^ above that 
of the best-fit model by these Icr noise errors. However, the the true errors on the size estimates may be 
somewhat larger due to the approximate modelling of the point-spread function and to asymmetry and 
substructure in the real radio galaxies. 

4.2 Results 

Firstly, the point-source models are generally a poor fit, giving a large x^- Using the estimated noise 

errors on x^, we can quantify the rejection of a point-source profile as ^ 2a for 6C1017+37, but ^ 6-20tT 
for the other 9 radio galaxies. Even for the QSO, a pure point-source is rejected by 6.7cr For all objects 
except 6C1256-I-36 and the QSO, there is a disk model and/or a bulge model with a x^ sufficiently low to 
indicate consistency within 2(t. In some cases, x^ is very similar for the best-fit disk and bulge models, 
for other galaxies one model may be favoured by as much as several cr. 

Figure 3 shows histograms of the area of each radiogalaxy above the 8 thresholds, together with 
the areas derived from the point-source model and the best-fitting disk and bulge models. To compare 
the models and data more directly, we also show radial intensity profiles. Figure 4 shows the mean 
if-band intensity in A(r) = 0.5 arcsec annuli about the centroid of each radio galaxy, together with the 
point-source profile and the best-fitting disk and bulge models. 

Table 2 lists the apparent surface brightness (SB) of each galaxy within the central 0.5 arcsec. 
This will correspond only approximately to the intrinsic properties of the galaxies, but does show some 
differences between them. The mean apparent central SB is 19.07 ± 0.13 K mag arcsec"^ for the 10 
6C radio galaxies, but the QSO is significantly higher in central SB and the galaxies 6C1017+37 and 
6C1129-I-37 are particulaly low. The morphologies of the individual galaxies are described below. 
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5. Morphology of Individual Radio Galaxies 

1. 6C0822+39 appears symmetric with a profile well- fitted by the bulge model, which is strongly ('^ Aa) 
favoured over a disk model. This galaxy appears in the if-band to be a normal elliptical. 

2. 600943+39 is small symmetric object, consistent with either disk or bulge at the la level. On the 
basis of this data, it could be another normal elliptical. 

3. 6C1011+36 is small and asymmetric with a much fainter secondary nucleus 1.7 arcsec from the main 

nucleus but clearly within the outer envelope, suggesting that it is currently merging with a much smaller 
galaxy. It may also be interacting with two close companions, which produce the peak in the profile at 
r ~ 4 arcsec. Neither a disk or a bulge profile is strongly favoured. 

4. 6C1017+37 is another small symmetric object, apparently the smallest of these galaxies. It is consistent 
with cither disk or bulge models at the Ic level, and the only galaxy consistent at the < 3(t level with 
being a point-source. The low central SB is presumably due to the object being faint and unresolved, 
rather than being an intrinsically low SB galaxy. There is no excess above the fitted models at either 
small or large radii, so this object could simply be a small elliptical. 

5. 601123+34 is round and symmetric, with a profile consistent with both the disk and bulge models. 
There is another galaxy about 5.5 arcsec away but no obvious indication that the two are interacting. 
This could be a normal elliptical. 

6. 601129+37 is extended and asymmetric with two intensity peaks of similar luminosity. These lie well 
within the same envelope, so that PISA detects 601129+37 as a single galaxy even with the deblend 
option. The outer regions appear disturbed with some evidence of trails, especially to the north. This 
is the only galaxy in this sample for which the profile significantly (~ 4(t) favours an exponential model 
over a bulge model. It also has a low central SB (0.49 mag less than the sample mean) which, as the 
galaxy is clearly resolved, would imply that the intrinsic central SB is lower than that of the other radio 
galaxies. A recent red-band WFP02 image of this galaxy, to be described in detail by Best et al. (in 
preparation), also showed that 

(i) the two peaks seen on the K image (centroids separated by 1.9 arcsec) lie within elongated 
structures with tails pointing inwards to the galaxy centre, suggesting that they are 'hotspots' associated 
with the outward passage of the radio jets rather than the two components of a merger. Olose to the 
centre the jets appear to be aligned with the 137 degree position angle of the two peaks of radio emission 
(AUington-Smith 1982; Eales et al. 1997), although the two radio lobes lie much further out, being 
separated by 14 arcsec. Further out from the centre, the ii-band jets curve slightly toward a W-E axis, 
causing the K-h&nd peaks to be slightly misaligned from the radio axis. 

(ii) The faint extension to the north is a small spiral galaxy of high i?-band (Arost — 3000A) SB, which 
appears to be interacting with the larger radio galaxy. AUington-Smith et al. (1982) had previously found 
601129+37 to be double on a ground-based i?-band OOD image, with the two components separated 
north-south so presumably corresponding to the radio galaxy and this companion rather than the two 
radio jet hotspots. 

(iii) No central point source is visible at Arest — 3000A. 

7. 601204+35 appears symmetric, slightly favouring (~ la level) a bulge over a disk profile. There is 
no excess above the fitted bulge model at large radii. This galaxy seems consistent with being a normal 

elliptical. 

8. 601217+36 is extended with a bulge profile favoured over an exponential disk by ~ 2a. The outer 
regions look slightly asymmetric, with some isophotal twist apparent on the contour plot. The radial 
profile is very close to the bulge model except for a small excess at r 5 arcsec due to a much smaller 
companion galaxy. This could, for example, be a normal elliptical just slightly perturbed by a near-miss 
encounter, resulting in a slight asymmetry. 

9. 601256+36 is extended and double with a less luminous secondary nucleus, 2.5 arcsec from the main 

nucleus. The outer regions are asymmetric and much more extended than in the other double galaxy 
(601011+36), with the second nucleus lying within what appears to be an inclined disk. The galaxy 
may also be interacing with its much smaller nearby companion (there a is hint of a bridge between the 
two). The galaxy is not fitted well by either a disk or bulge model, having too much area at the lowest 
thresholds. However, the radial profile looks closer to the bulge model at r < 2 arcsec, and the central SB 
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is typical of this sample (in contrast to 6C1129+37), suggesting that the primary galaxy is an elliptical. 
This is consistent with the visual impression of an elliptical merging with a large, lower SB disk galaxy. 

10. 6C1257+36 is closer to a bulge profile than a disk, although only by ~ la. The profile is well-fitted 
by the bulge model except for an excess at r 4 arcsec, due to a small companion galaxy. This galaxy 
could be another normal elliptical. 

11. The QSO 5C6.8 contains a very bright central nucleus, giving a central SB higher than that of any 
of the 10 radio galaxies, and a very small rhi- However, the object is not a pure point source. The 
K-hand image also shows the extended, lower surface brightness host galaxy, which appears disturbed, 
suggesting a recent merger or interaction. The profile is close to the fitted models at r < 3 arcsec but 
shows a significant excess at r ~ 4 arcsec, where it is more extended and resembles the profiles of the 
radio galaxies. Consequently, neither the disk or bulge models arc are a good fit, and a pure point-source 
model is also rejected. Two-component models and higher resolution data would be needed to investigate 
two-component systems of this type, by separately estimating a point-source contribution and an ru for 
the underlying galaxy (see Best et al. 1997). 

Despite the long wavelength of observation, these galaxies show a significant morphological diversity 
as well as a wide range of half-light radius. Two galaxies, 6C1011-I-36 and 6C1256-I-36 appear to be 
merging doubles. The 6C1129-I-37 galaxy is especially unusual - it is interacting and has luminous 
features associated with the radio emission visible on the ii'-band image. We refer to this for now as a 
'radio jet object'. 

Dunlop and Peacock (1993) and Ridgway and Stockton (1997) found that many z ^ 1 SCR galaxies 
do show an alignment between their K-hand isophotcs and radio axes, although no alignment was seen for 
Parkes radio galaxies an order of magnitude less radioluminous. Our 6C sample is intermediate in radio 
luminosity, so might be expected to contain at least one galaxy with visible if -band/radio alignment. 
The 6C1129-I-37 galaxy appears to be the only galaxy of these 10 with resolved aligned features, but 
it may also be significant that 6C0943-I-39, 6C1017-I-35 and 6C1204-I-35 show a close alignment (within 
10 degrees) of their JT-band position angles and radio axes (Eales et al. 1997). 

The 6C1129-I-37 galaxy is also the only object in this sample with a ii'-band profile which appears to 

significantly favour an exponential over a bulge model. Further observations may be needed to determine 
whether it is truly a disk galaxy or an elliptical with a profile distorted by the interaction and radio 
hotspots, but the low central SB compared to the other large double galaxy (6C1256-I-36), might favour 
a later Hubble type. 

The other seven galaxies appear at least consistent with the elliptical morphology seen for most radio 
galaxies, so will be referred to for now as 'bulge/generic' galaxies. We now consider the relation between 
angular size, luminosity and morphology for these and other radio galaxies. 

6. Size vs. Absolute Magnitude Relation 

We estimate the absolute magnitudes of the radio galaxies in the rest-frame Cousins i?-band (Acent — 
6560A), by subtracting the distance modulus from each observed K magnitude (Table 1), and adding a 
modelled correction to convert from the observer-frame Jf-band to the rest-frame i?-band. 

Roche et al. (1997) describe a model for the evolution of E/SO galaxies, in which star-formation begins 
at a high redshift, 16.5 Gyr ago, and decreases rapidly with an exponential timescale tsfr = 0.5 Gyr. 
As almost all star- formation occurs at z > 3, the evolution is passive at the redshifts we are concerned 
with here. The star-formation history is converted into an evolving spectral energy distribution using the 
stellar evolution models of Chariot and Bruzual (1993), and assuming a Salpeter IMF. Figure 5 shows 
the Rrest — Kobs correction, as a function of redshift, derived from this model's evolving spectral energy 
distribution, and also shows Rrest ~ -^o6s from much a bluer model in which star-formation begins at the 
same epoch but continues at a constant rate. 

At z — 1.1, the constant SFR model predicts much bluer observer- frame colours than the E/SO 
model, R — K = 3.5 compared to R — K = 5.5, but the Rrest — Kgbs correction differs by only 0.21 mag. 
At 2; ~ 1, many radio galaxies have R — K colours close to a passively evolving model, which essentially 
defines the red envelope of their colour-redshift distribution, but some of the more radioluminous (i.e. 
3CR) galaxies are as much as ~ 2 mag bluer and would be closer to the constant SFR model (Lilly and 
Longair 1984; Dunlop and Peacock 1993). Assuming an E/SO model should therefore give a conservative 
estimate of Mj^, accurate for the redder galaxies but an underestimate by as much as 0.21 mag for any 
objects which are as blue as the bluest 3CR galaxies. 
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Tabic 3 gives our estimates of Mr. All of these 6C galaxies are more luminous than the zero- 
redshift L* {Mji ~ —22.57 for ellipticals), with a mean M/? for the radio galaxies (excluding the QSO) of 
—24.25 ±0.18. To investigate the M^ — rhi relation, we adopt the disk model size estimate, r/i; = 1.68re,j.p 
for all objects with a smaller minimum for the disk model (6C1011+36, 6C1017+37, 6C1129+37, 
6C1256+36) and the bulge model estimate ru = Te for the others. Note that adopting instead the 
bulge-model rhi for the four objects with a lower for the disk-model would increase the errors on 
their rhi estimates, but would not significantly change their best-fit rui- If disk galaxy is observed at an 
inclined angle its apparent area will be reduced. To estimate the intrinsic size of inclined galaxies we 
would obviously use Thi on the major axis rather than the smaller Vhi derived from the apparent area. 
The observed area relative to that seen if the galaxy was seen face-on will be approximately (1 — e), where 
e is the image ellipticity as measured by PISA. We therefore correct the rui estimates for inclination by 
multiplying by (1 — e)~°-^. Table 3 gives the ellipticities and corrected rhi estimates, again with ±1(7 
errors, both in arcsec and in kiloparsecs for the assumed cosmology. The unweighted mean of the rni 
estimates for the radio galaxies (excluding the QSO) is 6.1 ± 1.3 kpc. 

Figure 6 shows rui against Mr for these galaxies, with symbols indicating morphological types, 

together with the 21 3CR galaxies from Best ct al. (1997), with the fitted bulge-model radii converted 
to kpc for our assumed cosmology. These 3CR galaxies cover a wide redshift range of 0.47 < z < 1.27, 
with Zmean = 0-89 Compared to our sample's Zmean = l-Hj so were divided into two subsamples, the 9 

galaxies at z > 0.93, with a mean redshift (z-mean = 1-10) alomst identical to our 6C sample, and the 
12 at 2 < 0.93, with 

Zmean — 0.73. The absolute magnitudes of the 3CR galaxies were estimated using 
the /T-band aperture magnitudes of Lilly and Longair (1984), extrapolated to total magnitudes assuming 
bulge profiles with r,. from Best et al. (1997), and adopting the same E/SO model Rrest — Kobs correction 
as used for the 6C galaxies. 

The SCR galaxies are on average more luminous (mean Mr = —24.88 ± 0.15) and much larger 
(mean rhi = 17.9 ± 1.5 kpc) than the 6C galaxies. The 3CR galaxies at 0.93 < z < 1.27 are more 
luminous (mean Mr = -25.30 ±0.16) than those at 0.47 <z< 0.93 (mean Mr = -24.56 ±0.19, but not 
significantly diff'erent in size (mean rhi = 17.5 ± 1.7 kpc compared to 18.2 ± 2.4 kpc). The 3CR galaxies 
at 0.93 < z < 1.27 are on average 1.05 ± 0.24 mag more luminous in the i?-band, and ~ 0.46 dex larger, 
than our 6C galaxies at the same mean redshift (note that this result will not depend significantly on the 
assumed Qq). This magnitude difference is greater than the 0.7 mag reported by Bales et al. (1997) for a 
z > 0.6 sample, probably as a result of our exclusion of the less evolved 3CR galaxies at 0.6 < z < 0.93 
and our use of total rather than metric magnitudes (giving a greater difference in luminosity between 
small and large galaxies). 

To compare these galaxies with those at lower redshift. Figure 6 shows local rhi — Mr relations, 

converted from the Vhi — Mb relations used by e.g. Roche et al. (1997) to the i?-band assuming the z = 
colour oi B — R ^ 1.44 for ellipticals (as given by the model) and B — R = 1.0 for spirals. 

For ellipticals, the Bingelli, Sandage and Tarenglii (1984) relation becomes 

log(rw/kpc) = -0.3(Mh + 20.01) 

for Mr < -21.44, and 

log(rw/kpc) = -0.1(Mfi+ 17.14) 

for Mr > -21.44. 

For spirals, the Freeman (1970) surface brightness becomes, 

log(rw/kpc) = -0.2Mij - 3.42 

Figure 6 also shows an observed sample of low redshift radio galaxies, from Owen and Laing (1989), 
who estimated absolute magnitudes in the Cousins R band and bulge-model-fitted radii for z < 0.2 
radio galaxies of a number of types. We plot their Classical Double FRII (Fanaroff and Riley 1974) 
radio galaxies, the type most appropriate for comparison with the higher redshift samples, with sizes and 
magnitudes converted to Hq ~ 50 km s~^Mpc~^ and their isophotal absolute magnitudes extrapolated 
to total magnitudes assuming bulge profiles with their fitted effective and isophotal radii. These 24 low 
redshift {zmean = 0.11) radio galaxies are less luminous (mean Mr = — 23.37±0.12) than those at higher 
redshift, and have a very wide range of sizes with a mean rhi of 12.5 ± 2.3 kpc, intermediate between the 
6C and 3CR galaxies. 
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The QSO is displaced far from tlie all other sources on this plot, on account of its central point- 
source. As we cannot measure the true galaxy size with this data, this object must be excluded from our 
discussion of rhi — Mr relations. 

The 6C galaxies are clearly shifted relative to the rhi — Mji relation of local ellipticals, in the direction 
of a higher intrinsic SB. As noted by Best et al. (1997), the SCR objects also tend to be enhanced in SB 
relative to local ellipticals. Note that estimates of the intrinsic SB will be independent of the assumed qo 
and Hq . The r^; — Mji relations of both the z ^ and z ^ 1 radio galaxies do not appear to follow the 
A(log rhi) = — 0.2A(M) slope of constant SB typical of spirals - they are noticeably steeper and more 
consistent with the A(log rhi) = — 0.3A(M) slope of the giant ellipticals. Hence to estimate the surface 
brightness evolution AR relative to local E/SO galaxies of the same size, we fit the function 

log(rM/kpc) = -0.3{Mr + AR + 20.01) 

For the 7 6C galaxies we classed as bulge/generic types, is minimized for AR = 1.79 ± 0.22 mag, with 
a of only 6.8 indicating that all 7 are consistent with this single rhi — Mr relation. One of double 
objects (6C1011+36) lies on the same relation. It appears to be undergoing just a minor merger, of a 
large elliptical with a much less massive dwarf galaxy, and this might not greatly affect the total size or 
luminosity. However, the other double object (6C1256+36) and the 'radio jet galaxy' 6C1 129+37 are of 
lower mean SB, corresponding to AR = l.Ol^to 13 mag and AR = 1.05^0 *17 m^-g respectively. The image 
of the 6C1256+36 galaxy was suggestive of an elliptical merging with a large, low surface brightness disk 
galaxy. This would obviously decrease the mean SB of the combined object and give it a more disk-like 
profile, while accounting for the central SB remaining high. In contrast, the radio jet object 6C1 129+37 
has a intrinsically low central SB as well as a low mean SB. 

An error-weighted fit to all 10 6C objects gives A_R = 1.26 ± 0.18 mag, but with a high of 62.16 
indicating that the galaxies are not all consistent with a single r;,; — Mr relation - there is a significant 
difference between the two relatively low mean SB objects and the relation defined by the other 8. 
An unweighted fit of the same relation to the 24 low redshift FRII galaxies gives AR = 0.12 ± 0.19 
mag, consistent with uncvolved ellipticals. An unweighted fit to the 21 3CR galaxies gives civolution of 
AR = 0.80 ± 0.15 mag, rather less than for our sample. However, if we separately consider the lower 
redshift (0.47 < z < 0.92) and higher redshift (0.92 < z < 1.27) SCR subsamples, the former gives 
AR — 0.51 ± 0.18 and the latter AR = 1.20 ± 0.18. There appears to be surface brightness evolution 
of 0.69 ± 0.25 mag between the respective mean redshifts of 0.73 and 1.10, with the higher redshift 
subsample being much closer in size-luminosity relation to the 6C galaxies at Zmean = l-H- We discuss 
the implications in Section 7 below. 

7. Discussion 

7.1 Surface Brightness Evolution 

Using ground-based data, we have carried out a preliminary investigation of the K-hand angular size 
and morphology of 6C radio galaxies at 1 < a; < 1.4. Our first result is that the 6C galaxies typically 
have smaller ii'-band half-light radii than SCR radio galaxies at similar redshifts. This appears to exclude 
the hypothesis that SCR and 6C sources at z ~ 1 have identical host galaxies and the higher K-hand 
luminosity of the SCR sources (e.g. Eales et al. 1997) is due solely to their more powerful active nuclei. 
If this was the case, the SCR galaxies would have smaller half-light radii than their 6C counterparts, due 
to the greater dominance of the central point-sources. The difference in K luminosity appears to be due 
primarily to the stellar components, which in the SCR galaxies are considerably larger and presumably 
more massive. 

We estimate that the Zmean = l-H 6C galaxies show a mean surface brightness enhancement of 
1.26 ± 0.18 mag in the rest-frame Cousins i?-band, relative to the rhi — Mr relation of local ellipticals. 
This degree of evolution is consistent with the 1.20 ± 0.18 mag wc estimate for a subsample of more 
radioluminous SCR galaxies at the same mean redshift. However, the 6C galaxies which most resembled 
normal ellipticals appeared to show greater SB evolution of 1.79 ± 0.22 mag. One possible explanation 
for this higher estimate is that some radio galaxies contain central point sources, which with the limited 
resolution of our data are not distinguished from the underlying galaxy but do cause an underestimate 
of rhi- This would cause the SB evolution to be overestimated, and if the effect on rhi is very large 
might even account for the difference in the mean rhi of the SCR and 6C galaxies. To estimate the likely 
effect of point-source contamination, we consider a bulge-profile model with rg = 1.8 arcsec - the mean 
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size of the Best et al. (1997) SCR galaxies - to which is added a central point-source contributing from 
0% to 100% of the total flux. These models were convolved with seeing profiles and normalized to have 
K ~ 17.4 with sky noise typical of the real Redeye data, and then analysed in the same way as the real 
radio galaxy images. 

Figure 7 shows the estimated bulge-model rg as a function of the point-source contribution to the total 
flux. A point-source component of 10% reduces the fitted by ^ 25%, so could increase the estimated 
AR by 0.42 mag. This would be sufficient to account for much of the difference in the estimated \DeltaR 
between the bulge/generic 6C galaxies and the 3CR galaxies. However, to reduce the fitted Vhi from 
1.8 arcsec to 0.55 arcsec, the mean size of the 6C galaxies, would require a much larger point-source 
contribution of > 30%. It is extremely unlikely that such luminous point-sources are typical of the 6C 
radio galaxies. Although we are at present only able to exclude a strong central point-source in the case 
of 6C1129-F37 (for which we have WFPC2 data), only 2/21 of the 3CR galaxies of Best et al. (1997) had 
profiles consistent with any point-source component of > 10% in the observed JC-band, and if anything 
the less radioluminous 6C galaxies would be expected to show less nuclear activity. 

Hence it seems that central point-sources would not make more than a small contribution to the 
difii'erence in rni between 6C and 3CR galaxies, but they might still be present in the bulge/generic 
galaxies at the < 10% level and be sufficient to cause an overestimate of the surface brightness evolution 
by up to ~ 0.42 mag. This would still leave at least ~ 1.37 mag of genuine evolution. This is still less 
than the surface brightness enhancement of ^ 3 galaxies on the Hubble Deep Field (e.g. Jones and 
Disney 1997; Roche et al. 1997b), so there is no doubt that such a degree of SB evolution is possible for 
large galaxies at z > 1. 

In constrast, the SB of FRII radio galaxies at Zmean =0.11 showed only 0.12 ± 0.19 mag of evolution 
relative to the local elliptical rhi — Mr relation, so these objects are consistent with being unevolved 
ellipticals or with the 0.12 mag of passive evolution predicted by our model at their mean redshift. This 
confirms that the assumed Vhi — Mr relation is appropriate for powerful radio galaxies at low redshifts 
as well as for normal ellipticals, and hence that the Vhi — Mr relation of FRII radio galaxies in particular 
does shift with increasing redshift. As would be expected, a Zmean = 0.73 subsample of 3CR galaxies 
showed intermediate evolution, estimated as 0.51 ± 0.18 mag. 

This progressive shift of the rhi — Mr relation with increasing redshift is similar to that seen for 
normal radioquiet ellipticals over a similar redshift range. Schadc ct al. (1997) found the rest-frame 
blue-band surface brightness evolution of a mixture of field and cluster ellipticals to be best-fitted by 
0.782 mag for qo — 0.5, or ~ 0.92^; mag for a low go- This is consistent with the passively evolving E/SO 
model described in Section 6, which predicts A(Mb) = —0.91 mag at z = 1. In the rest-frame Cousins 
red-band, the model predicts slightly less evolution of A(Mi{) = —0.71 mag, with 0.63 mag a,t z = 0.73, 
consistent with observations for the 3CR galaxies, and 0.75 mag at ^; = 1.11, a little less than observed. 
Hence there may be at least marginal indications that some z > 1 radio galaxies in both the 3CR and 6C 
catalogs undergo more optical evolution than expected from passive evolution. This might be explained 
as the more rapid fading of those components of the radio galaxy luminosity not produced by old stellar 
populations, e.g. young stars formed along the radio jets (see Best et al. 1996), and diffuse scattered 
light from the central AGN. Note also that a steeper IMF could increase the rate of passive luminosity 
evolution. 

7.2 Implications for Radio Galaxy Evolution 

The simplest evolutionary scenario is that surface brightness changes but does not. Although 
there is probably some increase with time in the mean size of galaxies, angular sizes from WFPC2 data 
suggest that, for ellipticals, this is a small effect out to z ^ 1 (e.g. Im ct al. 1996; Roche et al. 1997). 
Assuming no evolution in rhi, the estimated radii of the z ~ 1.1 6C galaxies would on the basis of the 
assumed local rhi — Mr relation indicate their local counterparts to have a mean absolute magnitude 
Mr = —22.6, i.e. about L* , and therefore suggest that many evolve into L ^ L* ellipticals. Most 
of the z < 0.2 FRII radio galaxies in the Owen and Laing (1989) sample are ellipticals of modest size 
and L ~ i* luminosity, although quite high radio luminosities of i(1400MHz) ~ 10^^-^ WHz~^. It 
is tempting to identify these as the low redshift counterparts of our 6C galaxies. Their typical radio 
luminosity corresponds to L(151MHz) ~ 10^''-^ WHz~^ for the average SED of luminous radio galaxies 
{a ~ —0.94), so this would require a factor of ~ 25 decrease in radio luminosity since z ~ 1.1. This is 
consistent with the 6C radio luminosity function undergoing the same r = 3.3 Gyr rate of evolution as 
estimated (Padovani and Urry 1992) for that of 3CR galaxies, but of course if our identification is correct 
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the evolution of the radio luminosity function must be due to the evolution of individual objects rather 
than simply a change in the number of galaxies which become powerful radio sources. 

Owen and Laing (1989) and Owen and White (1991) found that the classical double FRII radio 

galaxies at z < 0.2 were typically ^ 10^ '' higher in radio luminosity (Lrad) than FRI galaxies of similar 
optical luminosity, with the dividing line between the classes following a Lrad c< L^^^ relation (Ledlow 
and Owen 1996). As previously reported by e.g. Longair and Seldner (1979) and Prestage and Peacock 
(1988; 1989), the FRI radio galaxies tended to lie within clusters (on average Abcll class 0), while nearby 
FRIIs were usually in field environments. However, Yates, Miller and Peacock (1989) and Hill and Lilly 
(1991) found that powerful radio galaxies at higher redshifts {z ~ 0.5) were as likely to lie in Abell class 
0-1 clusters as in the field, suggesting a more rapid Lrad evolution for the sources in clusters. 

If any of these 6C galaxies lie in rich clusters, their low rcdshift counterparts might instead be FRI 
galaxies, like the smaller of those in the Owen and White (1991) sample, with much lower present-day 
radio luminosities of i(1400MHz) - lO^^-^ WHz'^ This corresponds to i(151MHz) ~ lO^''^-^ WHz"\ a 
factor of ^ 1000 down on the Lrad of our 6C galaxies. Yec and EUingson (1993) estimate that the optical 
luminosity evolution of AGN in rich cluster enviromnents (Abcll class ^ 1) is as rapid asT~1.0±0.2 
Gyr out to z ^ 0.6, but at higher redshifts must be similar to other AGN. If this form of evolution is 
paralleled by the evolution of Lrad it would predict a decrease in Lrad of a factor of ~ 1000 since z ~ 1.1, 
exactly the amount required. 

We estimate that the mean Mr of SCR and 6C galaxies at ~ 1.1 differs by 1.05 ± 0.24 mag. 
The mean radius of the SCR galaxies is much larger and with no size evolution would correspond to 
low-redshift counterparts with Mr = —24.2, several L* . On this basis we suggest that the difference in 
SCR and 6C luminosity at Xrest ^ is due to a difference in the stellar masses of the the host galaxies. 
As the radio luminosities differ by a factor ~ 6, this suggests a correlation Lrad oc M^^^n^^. 

A possible scenario explaining these observations is that; 

(i) At z > 1, the maximum radio luminosity is strongly correlated with the stellar mass of the host 
galaxy, approximately as Lrad oc M^. This relation could physically be related to the higher pressure 
of surrounding gas. and perhaps the more massive central black hole, in a more massive galaxy (see e.g. 
Bales 1992; Rawlings and Saunders 1991). At z > 1 many of the SCR and 6C galaxies may be observed 
while close to this upper limit, with the SCR galaxies typically being ~ 2.5 times more massive and hence 
~ 6 times more radioluminous. 

(ii) High redshift radio galaxies exist in a mixture of field and cluster environments. The peak Lrad 
attained during subsequent bursts of radio emission decreases more rapidly for radio galaxies in a rich 
cluster environment (Yee and EUingson 199S). By ^; 0, the cumulative effect of this environmental 
dependence amounts to a factor of ~ lO"*^'^ in Lrad between Abell class environments and the field, 
obscuring the Lrad correlation seen at z > 1. However, as cluster environment also influences the 
evolution of the radio morphology, a Lrad oc correlation reappears when nearby radio galaxies are 
separated into FRI and FRII types (Ledlow and Owen 1996). 

At any rate the high surface brightness and relatively small size of the 6C galaxies seems to argue 
that their evolution, at both optical and radio wavelengths, is no less rapid than that of the SCR galaxies. 
Although the K — z relation of 6C galaxies may have appeared consistent with no evolution (Eales et al. 
1997), their rhi — Mn relation clearly is not. Of course, the results of this paper, being based on a small 
dataset of marginally adequate depth and resolution, can only be regarded as a preliminary. We intend 
to investigate larger numbers of ^; ~ 1 radio galaxies, using NICMOS and ground-based adaptive optics 
to better constrain half-light radius and morphology, and to separate out bulge, disk, point-source and 
radio jet components. 

8. Summary of Main Conclusions 

(i) A sample of 6C radio galaxies at 1 < z < 1.4, observed in the /C-band, showed significant 
morpological diversity. Two appeared to be undergoing mergers, and another possessed two bright peaks 
aligned with the radio jets and may be closer to an exponential than a bulge profile, but the other seven 
were more consistent with normal elliptical morphologies in the i^-band. 

(ii) We estimated a mean half-light radius of 6.1 ± l.S kpc for the 6C galaxies, which is similar to 

that of present-day L ^ L* ellipticals and significantly smaller than the SCR galaxies studied by Best ct 
al. (1997). The difference between the measured sizes of 6C and SCR galaxies cannot be accounted for 
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by point-source contaniination, unless at the > 30% level, which seems implausible on the basis of the 
lack of point-sources in the 3CR galaxies. Hence the difference in Xrest ~ luminosity between 3CR 
and 6C galaxies at 2; ~ 1.1 must results from a difference in the sizes of the host galaxies, and not solely 
from a difference in the power of the central nuclei. 

(iii) The size- luminosity relation of the 6C galaxies at z ~ 1.1 is offset from that of either local 
ellipticals or FRII radiogalaxies at z < 0.2, indicating 1.0-1.8 magnitudes of surface brightness evolution 
in the rest-frame i?-band. This is similar to the 1.20 mag of surface brightness evolution seen for 3CR 
galaxies at the same redshift, suggesting that the 6C galaxies are undergoing a similar amount of evolution 
at optical wavelengths. 

(iv) The sizes of the 6C galaxies are consistent with their low-redshift counterparts being FRII radio 
galaxies like those studied by Owen and Laing (1989), a factor of ^ 25 lower in radio luminosity, suggesting 
radio luminosity evolution similar to the Padovani and Urry (1992) estimate for SCR galaxies. Any 6C 
galaxies which lie in clusters may undergo a much greater reduction in radio luminosity to become FRI 
objects. 
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Table Captions 

Table 1. The co-ordinates of the 11 radio galaxies, their if-band magnitudes as measured from this 
data, their spectroscopic redshifts, and the detection thresholds of the Redeye images. 

Table 2. A comparison of the areas of each radio galaxy above a series of 8 thresholds with those given 

by model profiles; (i) the for a point-source model; (ii) the and exponential scalelength (with ±lcr 
errors) of the best-fitting disk-profile model; (iii) the and effective radius (with itla errors) of the 
best-fitting bulge-profile model. 

Table 3. Estimated absolute magnitude in the rest-frame i^-band, the ellipticity of the radiogalaxy 
images as determined by PISA, and the half-light radius of the best-fitting model (with ±lc7 errors), 
corrected for inclination. 

Figure Captions 

Figure 1. Grey-scale plots of 15 x 15 arcsec areas of the -ftT'-band Redeye images centred on each 
radio galaxy North is at the top, east at the left. Top row (left to right): 6C0822-I-39, 6C0943-F39, 
6C1011+36. Second row: 6C1017+37, 6C1123-I-34, 6C1129-I-37. Third row: 6C1204-I-35, 601217-^36, 
6C1256-h36. Bottom row, 6C1257-h36, 5C6.8. 

Figure 2. Contour plots of the images on Figure 1. 

Figure 3. Histograms of the area (in 0.5 x 0.5 arcsec pixels) of each radio galaxy image above a series of 
8 intensity thresholds defined by PISA (the lowest being the detection threshold given in Table 1). The 
plots also show the areas measured above the same thresholds for a point-source (long-dashed) model, 
the best- fitting disk model (dotted), and the best- fitting bulge model (short dash), when convolved with 
the seeing profile and normalized to the same total intensity as the galaxies. 

Figure 4. Radial intensity profiles of the radio galaxies on the Redeye images, shown as the mean 
surface brightness in K mag arcsec"^ (with ±lcr errors) in annuli of 0.5 arcsec width around the image 
centroid, together with a point-source (long-dashed) model, the best-fitting disk model (dotted), and the 
best- fitting bulge model (short dash), all convolved with the seeing profile and normalized to the same 
total intensity as the galaxies. 

Figure 5. Modelled magnitude correction between the observer- frame JC-band and rest- frame Cousins 
i?-band, as a function of redshift, for a passively evolving E/SO galaxy model and a model with a constant 
star-formation rate. 

Figure 6. Estimated half-light radius (in kpc) against rest-frame i?-band absolute magnitude for the 6C 
radio galaxies, with symbols indicating morphology, together with the SCR radio galaxies from Best et 
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al. (1997), the low-rodshift FRII galaxies from Owen and Laing (1989), and assumed r^i — Mji relations 
for local ellipticals (dashed line) and spirals (dotted line). 

Figure 7. Bulge-model half-light radius (re) measured as described in Section 4 for a simulated Ve = 1.8 
arcsec bulge-profile galaxy to which is added an increasingly luminous central point-source. 
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